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Abstract. We investigate the big-bang nucleosynthesis in a Brans-Dicke model with a varying A term using the 
Monte-Carlo method and likelihood analysis. It is found that the cosmic expansion rate differs appreciably from 
that of the standard model. The produced abundances of 4 He, D, and barely Li are consistent with the observed 
ones within the uncertainties in nuclear reaction rates when the baryon to photon ratio r/ = (5.47 — 6.64) x 10~ 10 , 
which is in agreement with the value deduced from WMAP. 
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1. Introduction 

The standard model of the big-bang nucleosynthesis 
(SBBN) has succeeded in explaining the origin of light ele- 
ments 4 He, D, and 7 Li. Although the value of the baryon- 
to-photon ratio r\ has been derived from the observations 
of the Wilkinson Microwave Anisotropy Probe (WMAP) 
(Bennett et al. EZUU3"|l to be n 10 = 6.1±°* , the value 
seems to be inconsistent with the results of SBBN (Coc 
et al. I2UU4|I . Contrary to the excellent concordance with 
77 of WMAP for D, the abundance of 4 He by SBBN is 
rather low compared to that from WMAP. Therefore, non- 
standard models of BBN have been proposed with the 
Friedmann model modified (Steigman 2003). 

For non-standard models, scalar-tensor theories have 
been investigated (e.g., Bergmann 119681 Wagoner 119701 
Endo & Fukui IT3771 Fukui et al. IMTTjl . For a simple 
model with a scalar 0, it is shown that a Brans-Dicke 
(BD) generalization of gravity with torsion includes the 
low-energy limit string effective field theory (Hammond 
1996). Related to the cosmological constant problem, a 
Brans-Dicke model with a varying A(<^>) term (BDA) 
has been presented, and also investigated from the point 
of an inflation theory (Berman I1989|l . Moreover, it is 
found that the linearized gravity can be recovered in the 
Randall- Sundrum brane world (Garriga & Tanaka feOOOj) . 
Furthermore, scalar-tensor cosmology is constrained by % 2 
test for WMAP spectrum (Nagata et al. I2004[> where the 
present value of the coupling parameter loq = <jj(4>q) is 
bounded to be uj > 50 (4tr) and uj > 1000 (2a) in the 
limit to BD cosmology. 



In the mean time, BBN has been studied in BDA (Arai 

EMI 



et al. H5571 Etoh et al. HMTjl . A relation between BBN 
and scalar-tensor gravity is investigated with the inclusion 
of e + e~ annihilation to the equation of state, where the 
present value of the scalar coupling has been constrained 
(Damour & Pichon H999f) . On the other hand, it is sug- 
gested that a decaying A modifies the evolution of the 
scale factor and affects the temperature T r of the cosmic 
microwave background at redshift z < 10 4 , when the re- 
combination begins due to the decrease in T r (Kimura et 
al. l200ip . while a decaying A is found to be consistent with 
temperature observations of the cosmic microwave back- 
ground for z < 4 (Puy 2004). Therefore, it is worthwhile 
to check the validity of BDA related to the recent obser- 
vations. In the present paper, we investigate how extent 
BBN in the BDA model can be reconciled with r\ from 
WMAP. 



In §2, the formulation for BDA is given and the 
evolution of the universe in BDA is shown. Our re- 
sults of BBN are presented in §3 using the Monte-Carlo 
method (Cyburt et al. l2001[) . and constraints are given to 
the parameters inherent in BDA. We examine in §4 the 
evolution of the scale factor and the resulting abundances 
with taking into account the deviation from the equation 
of state p — p/3 during the stage of e + e~ annihilation. In 
§5, the likelihood analysis (Fields et AB. II 9961 is adopted 
to get most probable values and the accompanying errors. 
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Fig. 1. Evolution of the scale factor for BDA (// = 0.7, t = 13.7 Gyr, and rj 10 = 6.1) and Friedmann model. 



2. Brans-Dicke cosmology with a varying A term 

The field equations for BDA are written as follows (Arai 

et ai. rrwii : 



8tt„ lu 

1> 
1 



dcj) cj>- 



(1) 

(2) 



where lo is the coupling constant. 

The equation of motion is obtained with use of the 
Friedmann- Robertson- Walker metric: 



f dr 2 

ds 2 = -dt 2 + a{t) 2 I — 

l_ 1 — kr z 



+ r 2 d0 2 + r 2 sin 2 6 dcj) 2 \ , (3) 



where a(t) is the scale factor and k is the curvature con- 
stant. Let 1 be a scale factor normalized to its present 
value, i.e., x = a/ao, then we get from the (0,0) compo- 
nent in Eq. |QJ 



x\ k A uj 1 4>\ x 6 8ir p 

2 ---- I I I --- = — T> ( 4 ) 



x ) x 2 3 6 \ <j> I x cj> 3 



where p is the energy density. 

We assume the simplest case of the coupling between 
the scalar and matter fields: 



Y 2w + 3 p v ' 



(5) 



where fi is a constant. Assuming the perfect fluid for 
Eq. J2J reduces to 

dt^ x) = " 



r(p-3 P )x 3 



(6) 



2u + 3 

where p is the pressure. 

A particular solution of Eq. J2J is obtained from 
Eqs. (|TJ and ©: 



. 2?r — 1) _ 3 



(7) 



where /9 m o is the matter density at the present epoch. 
The gravitational "constant" G is expressed as follows 



2uj + 1 \ 1 
2u + 3^ ) 6' 



(8) 



The radiation density p r contains the contributions 
from photons, neutrinos, electrons and positrons at t < 
1 s. The total energy density is given as 



P 



Pm + Pr) Pr = Prad + Pv + Pe± ■ 



(9) 

Here the energy density of matter varies as p m — p m0 x~ 3 . 
The radiation density p r = p rCl x~ 4 except e + e~ epoch 
where e + e~ annihilation changes the relation T r ~ x^ 1 . 
We assume that the pressure satisfies p = p/3, which is 
legitimated only for relativistic particles. Then, Eq. © is 
integrated to give 



8717* \ 1 



(10) 



where B is a constant (Arai et al. H987[l . Although the re- 
lation p = p/3 does not hold during the epoch of e + e~ an- 
nihilation, as pointed out by Damour and Pichon (1999), 
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the inclusion of <f> measures small deviation from SBBN 
in the case of our interest, so that our solution can 
reasonably describe the evolution of 4> except the annihi- 
lation epoch. We consider that B affects the evolution of 
x significantly from the early epoch to the present com- 
pared to the contribution from e + e~ annihilation. As the 
consequence, neutron to proton ratio is affected seriously 
by the initial value of B. Considering the important con- 
tribution of p e ± to BBN, we examine the effects of e + e~ 
annihilation in §4. Hereafter we use the normalized values: 
B* = B/(10~ 24 g s cm~ 3 ), and r) W = 10 1( V The coupled 
equations Q), 10, and (|10J) can be solved numerically 
with the specified quantities: as for macroscopic quanti- 
ties, G = 6.672 x 10" 8 dyn cm 2 g" 2 , H = 71 km s^ 1 
Mpc" 1 (Bennett et al. l2UU5|> . and T r0 — 2.725 K (Mather 
et al. 11999(1 : as for microscopic quantities, the number of 
massless neutrino species is 3 and the half life of neu- 
trons is 885.7 s (Hagiwara et al. 2002). Though we adopt 
uj = 500, the epoch of the appreciable growth of \G/G\ 
is t < 10 3 s regardless of the value w (Arai et al. I1987[) . 
Therefore, even if we adopt a value to > 500 fWill I200T)| . 
we can get qualitatively the same conclusion by chang- 
ing the parameters /i and B* . We impose the condition 
|0/^| o = \G/G\ < 10" 13 yr^ 1 which is the most severe 
observational limit (Muller et al. I1991f) . 

We remark that BDA is an extension of the original 
form of BD and reduces to the Friedmann model when cj) 
= constant, fj, = 1, and We have A < if /x < 1, 

and (j>G > if both fi > 3 and Figure ^ shows the 

evolution of the scale factor for BDA with the relevant 
parameters in the present study and for the Friedmann 
model. Note that the difference in the expansion rate at 
t < 10 s in BDA. In particular, around t — 5 s, the curve 
x in BDA crosses that of the Friedmann model, which will 
have sensitive effects on BBN. Since A is proportional to 
p m o, p affects the evolution of the scale factor around the 
present epoch. In our BDA model, if \B*\ increases, the 
expansion rate increases at t < 10 — 100 s. It is remarked 
that the change in G between the recombination and the 
present epoch is less than 0.05 (2a) from WMAP (Nagata 
et al. I2004f) . which is consistent with BDA since |(G — 
G )/G | < 0.005 at t > 1 yr. 

3. Big-bang nucleosynthesis 

Changes in the expansion rate compared to the standard 
model affect the synthesis of light elements at the early 
era, because the neutron to proton ratio is sensitive to the 
expansion rate. 

For the BBN calculation, we use the reaction 
rates (Cyburt et al. I2U0T|> based on NACRE (Angulo et 
al. H999f) . We adopt the observed abundances of 4 He, D/H 
and 7 Li/H as follows; Y p = 0.2391 ± 0.0020 (Luridiana et 
al. 12003, D/H= 2. 78±g;f| x 10~ 5 (Kirkman et al. 127)173) . 
7 Li/H=(2.19 ± 0.28) x 10" 10 (Bonifacio et al.EJJJJ). 

Since the results of WMAP constrain cosmological pa- 
rameters, we calculate the abundance of 4 He, D and 7 Li 
paying attention to the value ry 10 =6.1. First, we carry 



0.3 r- 
0.28 - 




1 



Fig. 2. Light element abundances against r\ in BDA for 
/i = 0.7 and possible values of B*. 




Fig. 3. Same as Fig. [21 but for B* = and various values 
of \i. 



out the BBN calculations with use of the adopted experi- 
mental values of nuclear reaction rates given in NACRE. 
Figure H illustrates 4 He, D/H, and 7 Li/H for fi = 0.7. The 
abundance of 4 He is very sensitive to both B* and /i; it 
increases if \B* \ or /i increases. On the other hand, D and 
7 Li are more sensitive to /i than B* as seen from Fig. [3] As 
a result, 4 He and D/H are consistent with r/ obtained from 
WMAP in the range -0.5 < /i < 0.8 and -10 < B* < 10. 

Next, we perform the Monte-Carlo calculations to 
obtain the upper and lower limits to individual abun- 
dance using the uncertainties in the nuclear reaction rates 
(Cyburt et al. EOOTJ- Figure H illustrates 4 He, D/H and 
7 Li/H with 2cr uncertainties for B* = —2.5 and [i = 0.7. 
The light-shaded areas denote the regions of observed 
abundances, and the dark-shaded area indicates the limit 
obtained from WMAP. While the obtained values of 4 He 
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Fig. 4. Light-element abundances vs. 77 in BDA for /z = 
0.7 and B* = —2.5. Dashed lines show ±2<r uncertainties 
in nuclear reaction rates. The dark-shaded area indicates 
the constraint by WMAP and light-shaded areas denote 
regions of observational abundances. 

and D are consistent with 77 by WMAP, the lower limit in 
7 Li is barely consistent. 

4. Effects of e + e annihilation on BBN 

In the previous sections, we have assumed the equation of 
state p = p/3 in Eq. © to obtain Eq. (11 Oil at the epoch 
of e + e~ annihilation. Let us discuss the effects of e + e~ 
annihilation on the evolution of the scalar field and the 
scale factor due to the deviation from the relation p = 
p/3. The electron-positron pressure and energy density 
are written with the variable C = m e /k B T r as follows 



2m 



K 2 «) 



3p e 



2m e \ ^ 
^3 2^ 

n=l 



(11) 



(12) 



where ?i is Planck's constant in units of 27r, fc B is 
Boltzmann's constant, and m e is the electron rest mass. 
Ki (i = 1 and 2) are modified Bessel functions of order 
i (e.g. Damour & Pichon 1999). In the numerical calcu- 
lations, the summations in Eqs. l|llfl and (|12fl are taken 
over rt = 1 — 10. We can get the scale factor by integrating 
Eq. |TJ} with the aid of Eq. ©. To see the effects of e+e - , 
we take the form: 



8lTfl 

2w + 3 



(p e - 3p e )x 3 dt + B. 



(13) 



A direct comparison is made for the evolution of the 
scalar field. The results are shown in Fig. [3 where the solid 
line indicates the case of Eq. (|13|) with B* = —2.43 and 
[i = 0.7, and the broken line is the case of Eq. (|10|l with 
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Fig. 5. Evolution of the scalar field. The solid line refers 
the integration of Eq.(13) with B* = —2.43, and the bro- 
ken line is for Eq. (10) with B* = -2.50. 



B* = —2.50 and /j, — 0.7. These sets of parameters yield 
the same macroscopic quantities given in §2. Although we 
can appreciate the slight difference at t < 10 3 s, it remains 
small during and after the stage of BBN. The effects on 
the evolution of the scale factor are minor and the change 
in Y p is found to be at most 0.1% compared to that ob- 
tained in §3. We can conclude that since the effects of B 
in the range — 10 < B* < 10 is much larger compared to 
those of e + e~ annihilation, the deviation from the relation 
p = p/3 due to e + e~ does not change our results qualita- 
tively. However, we note that even small differences in Y p 
may affect the detailed statistical analysis combined with 
theoretical and observational uncertainties performed in 
the previous sections. 

5. Discussion and conclusions 

We have carried out the BBN calculations in the fj, — B* 
plane and obtain the ranges —0.5 < [i < 0.8 and —10 < 
B* < 10 that are consistent with both the abundance 
observations and 77 obtained from WMAP. 

To evaluate uncertainties of theory and observations, 
we calculate normalized likelihood distributions in BBN 
(Fields et a UTMHl Hashimoto et al. 12003 • In Fig. we 
show the likelihood functions for 4 He, D, and 7 Li. The 
combined distributions, L A7 = i 4 • L 7 and £ 2 47 = ^2 ' 
i 4 • L 7 are shown in Fig. [7J As the result, we get the 95% 
confidence limit of m 5.47 < 77 10 < 6.64. 

The consistency holds within la error for He and D, 
and 2cr for 4 He, D, and 7 Li. Though new reaction rates 
recently published (Descouvemont et al. l2004|l will change 
the errors to some extent in the likelihood analysis, our 
conclusion holds qualitatively. 

Our previous studies (Etoh et al. 119971) showed 1 < 
/i < 3 if A > for large value of u>. In the present case, 
the A term becomes negative from Eq. (JTJl for fi < 1: this 
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Fig. 6. Likelihood function as a function of r)io for 
4 He (i 4 ), D (L 2 ) and 7 Li (L 7 ). The vertical lines indi- 
cate upper and lower limit to i] by WMAP. 




Fig. 7. Combined Likelihood function for two (L 47 ) and 
three-elements (L247) • 



would not conflict with available observations and/or ba- 
sic theory (Vilenkin 2004). Alternatively, if we consider 
A = A + A((/>) with |A(0 o )/A o | < 0.01, then the cosmo- 
logical term becomes consistent with the present observa- 
tions. Though the evolutionary path in the early universe 
can deviate from the Friedmann model (Arai et al. [1987) , 
parameters in BDA must be searched in detail for val- 
ues of u> > 500 to get quantitative results of BBN. We 
note that it is shown that negative energies are present in 
scalar-tensor theories, though it is not clear how to iden- 
tify them definitely (Faraoni[2P04). 

To save the apparent inconsistency for SBBN, effects of 
neutrino degeneracy, changes in neutrino species, or other 
new physical processes have been included (Steigman 
2003 ) . In our model, we need only a scalar field that could 
be related to a string theory fHammond ll996|l . It is noted 
that the original BD cosmology (fi = 1) would be limited 



severely by the more accurate observation of light elements 
and/or the future constraints for rj as shown in the present 
investigation. 
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